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A new functional for the entropy that is asymptotically correct both in the high and low density 
limits is proposed. The new form is 

S = S {ld) + S {ln) + S {r) + S (c) 

where the new term 8^"' depends on the p-bodies density fluctuations a p and has the form 

S (c) = <iV} J In 2 - 1 + V ^-^a p - [exp(a 2 - 1) - a 2 ] I + S 

I U pl J 

S renormalizes the ring approximation S^ r '. This result is obtained by analyzing the functional 
dependence of the most general expression of the entropy: Two main results for S^"' are proven: i) 
In the thermodynamic limit, only the functional dependence on the one body distribution function 
survives and ii) by summing to infinite order the leading contributions in the density a new numerical 
expression for the entropy is proposed with a new renormalized ring approximation included. The 
relationship of these results to the incompressible approximation to entropy is discussed. 
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I. INTRODUCTION 



Entropy is one of the very important and challenging thermodynamic quantities in statistical mechanics; this is so 
because it depends on all the n-particlc distribution functions. The problem is not the lack of exact expressions but to 
derive equations that are both accurate and manageable from a theoretical and numerical point of view. Amongst the 
exact expressions we cite the classical work of Nettleton and Green 13], and, more recently, while approximate 

expressions can be found in || ^] . In Q we have shown how an infinite subset of terms (dependent exclusively on the 
one- and two- body distribution functions) can be analytically summed giving rise to the so-called ring approximation 
and, through a minimization of the free energy functional, the well known HNC approximation la] is obtained as 
an optimized superposition approximation. Later on, Bush et al || derived sets of integral equations by analyzing 
several levels of approximation to the grand potential function and, in a recent and very interesting article, Puoskari 
|p| extends the ring approximation to three particle functions showing how the HNC2 equations either of the Wertheim 
|)| or the Baxter Jicfl variety can be obtained. Baranyai and Evans || showed that, even though the derivations are 
done in the grand canonical ensemble, the entropy equations are, in fact, ensemble invariants if local expressions are 
used for the entropy. In this way, the comparison with canonical ensemble numerical simulations is justified. They also 
analyze the convergence range of some needed integrals showing how this range increases at high densities. Wallace 
also worked with the same type of expression and, by analyzing the behaviour of density fluctuations, proposed an 
incompressible approximation for the entropy of dense fluids (which, as Wallace himself states, is decidedly wrong in 
the low density regime) . He then concluded that the dilute gas and dense liquid regimes occupy different regions of the 
phase space. These conclusions have been criticized |J . Laird and Haymct |p"l|| have extended the ring approximation 
to mixtures and applied it to electrolytes. They found that the correct Debye-Huckel expression for the entropy in 
the low concentration limit is obtained when the ring approximation is included. They have also discussed |12| the 
incompressible approximation in dense fluids, proposed an expression that differs from the one of Wallace, and applied 
it with good numerical results. 

Summarizing, these works show that, at low densities the inclusion of the ring approximation gives a very accurate 
entropy equation and the incompressible approximation is hopelessly wrong, while, at high densities, the ring term 
overestimates the entropy and the incompressible approximation is reasonably accurate. It has been suggested that 
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more than two bodies correlations must be incorporated in order to have an accurate expansion and, in this respect, 
Puoskari's work || is quite promising. It is the purpose of this article to elaborate on the compressibility related 
contribution to entropy and improve the entropy expansion when truncated to the pair distribution level. In Section II 
we discuss the conceptual structure of the entropy when written as a functional of the n-bodies distribution function. 
In Section III it is shown by functional differentiation that the compressibility related contribution only depends on 
the one body distribution function as well as on thermodynamic parameters and, in Section IV, by summing three 
subsets of terms, we shed some light on the nature of the incompressible approximation and show how the correct 
low density limit can also be obtained. In Section V we present our conclusions and propose a new entropy equation. 



II. ENTROPY STRUCTURE 

We can think on two different criteria for the analysis of the entropy. The first one can be called the functional 
criterium; is the one we take when interested in the entropy dependence on distribution functions, e.g., when a 
variational principle is formulated. The second point of view is the numerical one, i.e. the main goal is numerical 
accuracy. It is clear that it is not necessarily true that the same expression will fulfill satisfactorily both goals. In 
particular, the incompressible approximation (see eq. ([h] )) is numerically correct at high densities, but has the wrong 
behaviour at low densities and also lacks a sound theoretical foundation. 

When the entropy dependence on all the n-bodies distribution functions is explicitly written, we obtain an approx- 
imate expression with the following distinct structure contributions : 
i) The ideal gas contribution 



-r- = W 



(1) 



where p is the number density and A the thermal wavelength, 
ii) The ever present logarithmic contribution 



S Qn) 



= -E/ rf{pK({p}H({p» (2) 

P >i J 



where n p ({p}) is the p-particles distribution function and w p ({p}) the irreducible p-bodies contribution to the po- 
tential of average force. More specifically, the link with the more usual notation is 

e ^ p ({p}) = n P ri{p-2}c{p} n P~2 Il{p-4}c{p} n P~4 ■ • ■ ,^ 

ri{p-i}c{p} n p-i ri{p-3}c{p} n p-3 ■ ■ ■ 

As usual, we have that g% = 1 + h% and, through the use of the generalized superposition approximation (GSA) 
we can write that, e.g. 

.g 3 ({3}) = [l + A 3 ({3})] [] .9 2 ({ 2 }) 

{2}C{3} 

, r , n M . A ,,.,», ri{3}C{4}53({3}) 

94 {4} = [1 + A 4 ({4})] 1 = (5) 

11{2}C{4} 52({2}) 

[1 + A 4 ({4})] J] [1 + A 3 ({3»] J] g 2 ({2}) 

{3}C{4} {2}C{4} 

which introduces the family of A p functions that, when different from zero, correct for the difference with the GSA. 
They can also be written as 

A p ({p}) = e^«P» - 1 (6) 

This S^ n) contribution is the one that, when functionally differentiated, gives rise to the In g contribution in the 
integral equations 
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iii) The ring term, which in its simplest, two-body version, is 

_ 1 ^2 

p>3 



IT^E^— J d {P}f[Mi)h 2 (12)h 2 (23)...h 2 (pl) (7) 



and it can be summed in homogeneous systems Q . This term is the responsible of the contribution h 2 — c 2 in the 
integral equations. The three bodies version is derived in B. 
iv) The compressibility related contribution 



k 



El/ d{p}f[ ni (i)A p ({p})r p ({n p }) 

p> 2 P' i=l 



(8) 



r ( < X \ - n{ P ~i}c{p} 9 P -i ri{p-3}c{p} 9 P -3 ■ ■ 

1 P\V L PS) — T-r T-r l y J 

ll{p-2}c{p} 9p~2 ll{p-4}c{p} 9p-i ■ ■ ■ 

Its first term is essentially the compressibility (A2 = h 2 ) and the sequence of products in eq. (0) stops when reaching 
either 173 or g 2 . As far as we know there are no previous studies of the whole series given in eq. (Q); the compressibility 
approximation focuses on the first term of this series, which for a one component homogeneous system is 

^ = Wf/ dvh 2 {r) = ^ (-1 + 02) (10) 
As in the dense liquid limit is a 2 <C 1, the incompressible approximation considers a 2 = in the whole p — T space. 

III. FUNCTIONAL DEPENDENCE 

Here we prove that, in the thermodynamic limit, all the functional derivatives of the compressibility contribution 
with respect to the distribution functions can be summarized in the equation 

'" S ""' // '' - Si P + 0(e- <N> ) (11) 



Sn p ({p}) 
Therefore, 

= " / d{l}ni({l}) + C(p, T) + 0(e-< N >) (12) 

C is an integration constant as far as the functional integration refers but, in fact, it depends on p 1 T. 

The derivation is straigthforward. The origin of the compressibility term is quite clear and eqs. (38-41) of Ref. 
H are the equations to look at. Equation (38) is our eq. (]l(]) and in eqs. (39-41) we see that each one of them has, 
amongst other terms, the integral j Y[ n i9pd{p}- When the GSA for g p is used (eq. the integral decomposes 

into a sum of two integrals J Y[ n i^p [1 + A p ] d{p}. The term without A p contributes to the ring approximation plus 
neglected terms (such as those shown in Ref. j|) and the term with A p contributes to S^- c \ 



As T p can also be written as 



r ( < X \ - n{p-i}c{ P } n v~i ri{p-3}c{p} n p-3 ■ ■ 

ll{p-2}c{p} n P-2 ll{ p _4}c{p} n P-4 ■ • 

using eqs. (§) and (|) we conclude that 



p>2 * 



ri{p-i}c{p} n p-i ri{p-3}c{p} n p~3 ■ ■ ■ 

Il{p-2}c{p} n P~2 ]l{p-4}c{p} n p-4 ■ ■ ■ 



(13) 



It is somewhat clear that each one of these integrals is related to p-bodies density fluctuations but a clearer explanation 
is to be found in the next section. This explains the origin of naming this contribution as compressibility related. 
Written in this way it is straightforward to show that the functional derivatives are 
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SS^/k 
<5ni(x) 



p>i 



•0(e 



-<AT) 



(14) 



SS^/k 
<5n s ({x s }) 



p>0 



■ W 1 



1+0 (m) 



O (e- <">) 



(15) 



and we arrive to eq. (Jll 

This result shows that, in the thermodynamic limit, the compressibility term does not contribute to any set of 
equations we may derive by functional diferentiation of a functional that includes the entropy; it only contributes to 
the constraint of fixed density. Therefore, if we are after a set of equations which are the consequence of a variational 
principle, we can rightly put the compressibility term aside from all the others, while if we are after a numerical 
approximation to the entropy we can assume, on physical grounds, absolute convergence of the whole entropy series 
and feel free to mix the compressibility term with all the others if numerical convergence is improved. 



IV. SERIES SUMMATION 

For the sake of this numerical goal we will cut the GSA to third order; in this way g p can be written in two equivalent 
forms 



n (M\ J 1 + £{2}C {P } MW) + £ { 3}C {P} MW) + . . . h p ({ P }) 

■ MiPlj _ i n {3} c {P} [i + a 3 ({3»] n {2} c {P } t 1 + Mm)] 



(16) 



In eq. (|l^) for the compressibility contribution we will sum to infinite order three subsets of terms. These subsets 
are clearly identified in the p — 3 summand of eq. (|l3|), i.e. 



S. 



(c) 



3 

| J d{3>nni(i) 
i J d{3}f[ ni (i) 



5s({3})- [] [! + ^({2}) 

{2}C{3} 



fo({3}) - II ^( ik ) - ^(12)^(13)^(23) 

i—1 k^i 



(17) 



A) The first subset includes the contribution of the integrals J h p d{p},p > 2. The series is 

c( c ) 00 i r p °° in \ 

^=Ei/^p}n^i)M{p})=E^ 

p=2 F ' J i=l p=2 F ' 

The moment-cumulant relation Q is 

M M k 

C M ({M}) = h M ({M}) J] m(i) = Y, {fc{m,}M} (-l) fc_1 (* - I)' II n ™* 



(18) 



(19) 



i=l 



k=l 



i=l 



Here, the partition of the coordinate set {M} in k disjoint subsets {m^M, 1 < « < fc is symbolized by {fc{rxii}M} and 
therefore X)fc=i {^{ m i}M} indicates the sum over all the partitions in k subsets and for each k is 1 < i < k. In this 
way (C p ) is related to the integrals {n Pi ) = J dpi"a Pi . On the other hand, (n Pi ) can be expanded in terms of (iV fc ) jl4|] 



(n p ) = {N(N -l)...(N-p+l)) = J2 «(P, k) (N k ) 



(20) 



k=l 



for p > 1, where s(p,k) are the Stirling numbers of first kind. One of its definitions is that (— l) p k s(p,k) is the 
number of permutations of p elements which contain exactly k cycles. They satisfy the recurrence relation 

s(p + 1, k) = s(p, k — 1) — ps(p, k), 1 < k < p 
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with starting values 

s(p,0) = 8(0, k) = S „ 
We also define the r-bodies density fluctuations a r by 

f 1 r = 1 / X 

Or = < ((N-<N» r ) , (21) 

I <7V> ' 1 

The first few (C p ) are then expressed in terms of the a r and Stirling numbers as 

2 

(C 2 ) = (N) (-1 + a 2 ) = (AT) ^ s(2, (22) 

3 

(C 3 ) - (JV) (2 - 3a 2 + a 3 ) = (iV) ^ s(3, i)oi (23) 

i=i 

4 

(C 4 ) = (JV) (-6 + lla 2 - 6a 3 + a 4 ) + 3 (TV) 2 a 2 , = (JV) ^ s ( 4 ; *)«< + 3 ( N f a l ( 24 ) 

i=i 

The 3 (N) 2 a 2 , term and similar ones from higher order (C p ) will be included in the next partial sum. In order to sum 
to infinite order the contribution of each r-bodies density fluctuations we need the result 



Therefore, the a r contribution is 



fe! t! 



rQ - - ^ (k?2) p /r! r > 2 (25) 



and our first partial sum, which includes contributions to infinite order of all the a r , is 

f-§^- W {^-i + §fif^} (-) 

Let us remark on some characteristics of this result: i) In 2 — 1 is the contribution in the absence of density 
fluctuations and it is also its high density limit (a p <C 1), ii) as ol p — ► 1 when p — > 0, eq. (|2q ) vanishes in the low 
density limit, iii) the series is rapidly convergent. Therefore, this sum goes in the right direction to improve on the 
incompressible approximation, both in its numerical results and in its theoretical foundation. This analysis makes also 
clear why this contribution is referred to as compressibility related. Lastly, the a p are easily expressed as integrals of 
the correlation functions; the first ones are 



«2 = l + ^y| d{2}]Jn 1 (i)h 2 ({2}) 

a 3 = -2 + 3a 2 + -L J d {3} JJ m(i)ft 3 ({3}) 



B) When in the p-th term of eq. ( fL3[ ) the expansions given in eqs. (jig ) are inserted, each hk, k < p expanded 
in terms of h 2 , A 3 and terms like the 3 (N) a 2 , which were left aside in the first partial series included, then all the 
unconnected (in the graph theory sense) terms cancel out and the first two sets of connected diagrams are those 
depicted in eq. (|l7|) . We first evaluate the sum of " star " products of h 2 bonds 

o(c) oo oo , p p „ 

^ = E^ = -E-E/ d{i}m(i)l[ d{k} ni (k)h 2 (ik) (27) 

p=3 p=3 i=l J kj=i J 

Each summand is easily evaluated as 
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and the second partial sum is 



g(.o) 

-L- = -(N) [exp(a 2 - 1) - a 2 ] (28) 



Its low and high density limits are and — e _1 respectively. 

C) This series is a sum of rings very similar to eq. (Q), its first term is given in eq. (17) and, as the symmetry 
number of p-rings is 2p, it can be written as 

^r = -\Y.\ /rf{p}f[ni (1)^(12)^(23)... fe(pl) (29) 

which can be summed for homogeneous systems in the same way as the original ring approximation was Q 

2 



<? (c) (N) r \ / \ (ph 2 (k) 

\ = / dk { In (1 - ph 2 (k)) - ph 2 (k) - V - ' ) CUD 

where h 2 (k) = J dvh{r) exp(27rzk.r) is the Fourier transform of h 2 (r) and the integration is over the k— space. This 
contribution can be added to the original ring approximation giving a renormalized ring approximation , which is 

- - (31, 



k 2 . 



5^- / i{2p}JJm(i)/i a (12)^2(23)...fta(2pl) 

P >2 P J i=l 



a sum over all even order rings and, for homogeneous systems, the result is 

~~k~ ~ ~p 



^ = M / dk/ln (l - (^ 2 (fc)) 2 ) + U 2 (fc)) 2 [» (82) 



V. CONCLUSIONS 



Collecting together the different results obtained, i.e. eqs. (p6|), (|28|) and (|32j ) with eqs. (|1|) and (|2[) we arrive to a 
new entropy equation which includes a partial summation of the compressibility related contribution, i.e. 



^-ln(pA 3 ) 



^ / d{p}n p ({p})w p ({p}) 
^ J dk{\n{l-(ph 2 (k)f)+(ph 2 (k)) 2 } + 



(33) 



(JV> ln2-l + £ 



(ln2f 



p=2 



Op? - (AT) [exp(a 2 - 1) - a 2 ] 



If, for just a moment, we neglect the renormalization in the ring approximation, the contribution due to the r-bodies 
density fluctuations gives not only a theoretical understanding on the nature of the incompressible approximation, 
but also a description that is essentially correct in the low and high density limits (0 and In 2 — 1 — e _1 respectively). 
As the ring approximation grows quite steeply when the density increases Jll],[l2| , its renormalization should have the 
right asymptotic behaviour. Let us also mention that these results extends trivially to mixtures (see, e.g. fl^l ) and, 
in this case, it is more convenient to work with the entropy per unite volume. 

We have also shown that, in the thermodynamic limit, the compressibility term only depends (as a functional) on 
the one body distribution function. Therefore, this functional dependence is such that it only enters in the constant 
density constraint and, in this way, the conceptual structure of the equation for the entropy is significantly simplified. 
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As our results apply to the full entropy functional, they are valuable to any functional minimization as, e.g., those in 
H,[| . It can also be mentioned that this theorem does not conflict with Laird-Haymet . They obtained the correct 
Debye-Huckel low density expression for the entropy by including the S^ d \ S^ ln \ S( r > plus the compressibility related 
contribution of eq. |Io| ). As this term and, in fact, all the sums we did, vanish when p — > (including the one that 
renormalizes the ring approximation), there is no contradiction between ours and theirs results. Lastly, as this result 
does not depend on the potential, it is also valid for the associative Wertheim-Ornstein-Zernike equation Jll|[l6| as 
well as systems with directional forces fll7|] . 

This new functional provides a robust and systematic way to develop fully analytical theories of liquids [fl8|-po[ , 
which will be examined in future work. 
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